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Abstract 

We report the growth and characterization of single crystals of iron chalcogenide superconductor 
FeCro.o2Se. There is an enhancement of the superconducting transition temperature (T c ) as com- 
pared to the single crystal of the parent compound Fei+^Se by about 25%. The superconducting 
parameters such as the critical fields, coherence length, penetration depth and the Ginzburg-Landau 
parameter have been estimated for these single crystals and compared with related Fe-based super- 
conductors. Analysis of the critical current data suggests a 51 pinning mechanism in this material. 
Thermally activated transport across the superconducting transition in the presence of external 
magnetic field suggest a crossover from a single vortex pinning regime at low fields to a collective 
flux creep regime at higher magnetic fields. Nature of charge carriers in the normal state estimated 
from Hall effect and thermal transport measurements could provide crucial information on the 
mechanism of superconductivity in Fe-based materials. 

PACS numbers: 74.25.Dw, 74.70.Xa 



3 

CZ2 



* Corresponding Author 

Email: ajay.thakur@iitp.ac.in 
' Corresponding Author 

Email: tomy@phy.iitb.ac.in 

1 



Iron chalcogenide a— FeSe crystallizes in the hexagonal NiAs type crystal structure with 
PQz/mmc space group, and is non-superconducting. A slight excess of Fe stabilizes the 
superconducting tetragonal phase (3— Fei +x Se (anti-PbO type; space group PA/nmm) with 
a superconducting transition temperature (T c ) of ~8.5K [1]. FeSe and its derivatives (com- 
monly known as the '11' system among the Iron based superconductors) have drawn im- 



mense attention 



as 



akin to the remarkable similarity in their Fermi surface with that of 
the FeAs based iron pnictide superconductors j^]. Therefore, despite their structural and 
compositional simplicity compared to the pnictide counterparts, they are conceived to be 
promising model systems to understand the physics of Fe-based superconductors. Struc- 
turally, FeSe comprises of stacked layers of corner sharing FeSe4 tetrahedra similar to the 
FeAs based materials. However, the spacer layers are absent. Analogous to the members of 
the pnictide family, there is a distortion in the FeSe layers across the structural transition. 
The properties of off-stoichiometric (3— Fei +x Se (or, (3— FeSei_. x ) was observed to be very 
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sensitive to the ambient pressure. Whereas, Mizuguchi et al. [5] reported a T c of 27 K at 

1.48 GPa in (3— Fei +x Se, Medvedev et al. [6] found a non-monotonicity in the evolution of 

T c with pressure with a maximum of 36 K at 8.9 GPa. The T c was seen to decrease within 

the pressure range 8.9-32 GPa. A non-monotonic variation of T c was also observed in the 

high pressure studies by Margadonna et al. |7J where T c was seen to peak at 37 K at 7 GPa. 
rn 

Garbarino et al. [8] on the other hand showed a monotonic increase in T c with pressure, 
where, an orthorhombic high pressure phase was seen to develop at 12 GPa with a T c of 
34 K at 22 GPa. Although the precise nature of the superconducting phase and the exact 
role played by external pressure is not very clear, the sensitive dependence of T c on ambi- 
ent pressure suggests a possibility of increasing T c via the chemical pressure route. Several 
groups around the world have tried different substitutions at both the Se site 0, [lo| and the 



Fe site 
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14J however, there have been no reports suggesting an enhancement of T c . Upon 



Ni substitution at the Fe site, the superconducting volume fraction was seen to enhance, 
but there was no concomitant increase in T c 12] . We have recently reported a route to en- 



hancement of 
samples 
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by the substitution of Cr instead of excess iron in polycrystalline FeCr^Se 
-Il8|. We have now been successful in growing single crystals of FeCr .o2Se. In 
this paper, we report the detailed superconducting and transport properties in these single 
crystalline samples. 

The single crystals used in this work were prepared using the self flux method for crystal 
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growth. Powders of 4N Fe, Se and Cr obtained from Alfa Aesar were homogenized in an 
agate mortar in the required stoichiometric composition (FeCro.o2Se), sealed in an evacu- 
ated quartz tube (10 -6 mbar) and preheated at 1050°C for 24 h. It was cooled down to 
800°C at a rate of 2°C/h and then the furnace was switched off. The sample thus obtained 
was annealed at 360°C for 36 h followed by quenching in liquid N 2 . The X-ray diffraction 
(XRD) and energy dispersive X-ray analysis (EDXA) were performed to confirm the struc- 
ture and elemental composition. Magnetization measurements were performed using the 
SQUID- Vibrating Sample Magnetometer (SVSM), Quantum Design Inc. (U.S.A.). Elec- 
trical and thermal transports were performed using the Physical Property Measurement 
System (PPMS), Quantum Design Inc. (U.S.A.). The scan amplitude was chosen to be 
2.0 mm and 1.0 mm for the PPMS-VSM and SQUID- VSM, respectively. 

The single crystals of FeCr .o2Se reported in this paper have layered planes held to- 
gether by weak van der Waals interaction and can therefore be cleaved easily similar to the 



FeSei-aTes UJ, NbSe 2 
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22j and Bi-2212 
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24] single crystals. The X-ray diffraction 



pattern obtained for one such cleaved single crystal piece is shown in the inset panel of 
Fig. 1 (a). Prominent (101) reflections are observed indicating that the (lOl)-axis of the 
single crystal is perpendicular to the cleaved surface. The panel (a) of Fig. 1 shows the tem- 
perature variation of the real part of AC susceptibility data. The onset temperature for the 
diamagnetic transition is marked as T° n; dta . In panel (b) of Fig. 1, we present a comparison 
of the DC magnetization and the electrical transport measurements across the supercon- 
ducting transition in FeCro.o2Se. Also marked in Fig. 1 (b) are the onset (T° n = 14.8 if), 
mid-point (T c "" d = 12.2 if), offset temperatures (T£** = 11.1 if) and zero resistance tem- 
perature (T c ° = 10.5 if) which are respectively the temperatures at which p reaches 90%, 
50%, 10% and 0% of the normal state value (pn) measured in the absence of an external 
magnetic field (H = 0). The temperatures for the onset of diamagnetism (T° n; dia ) and zero 
resistance (T c °) coincide with each other (hallmark of superconductivity). This temperature 
is identified as the superconducting transition temperature (T c ) and has a value of 10.5 K 
for our crystals of FeCr .o2Se. 

Any practical application of a typical type II superconductor demands an understanding 
of its vortex phase diagram [25] . In particular, it is of utmost importance to obtain the lower 
critical field, H cl and the upper critical field, H c2 values at various temperatures below T c . 
In panels (a) and (b) of Fig. 2 we present the H C 2 data obtained from the electrical transport 
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data via the magnetic field dependence of the T° n , T™ ld and values (see Fig. 1(b)) for 
both H || (101) and H _L (101). The value of the slope (^§fr-) are extracted from graph as -23 
kOe/K, -26 kOe/K and -25 kOe/K, for data corresponding to T c on , T c w and T c o// (identical 
within the first decimal place for both H || (101) and H _L (101)), respectively. The upper 
critical field values at zero temperatures H C 2(0) are therefore estimated as 24 kOe, 22 kOe 
and 17.5 kOe (for both H || (101) and H _L (101)) using Werthamer-Helfand-Hohenberg 
(WHH) formalism [26], where, H c2 (0) = -0.693 {^f)T c for the three different criteria, 
respectively. This points to an absence of magnetocrystalline anisotropy in single crystals of 
FeCro.o2>Se. The estimated values for H C 2(0) are very close to the Pauli paramagnetic limit 
Hp(0) = 1.84T C {27)], suggesting spin-paramagnetic effect as the dominant pair preaking 
mechanism similar to the claims made for the FeSe system with S and Te substituted at the 
Se site 0, Q. An estimate of the superconducting coherence length £(0) was made using 
the Ginzburg-Landau expression H c2 (0) = $o/27r£ 2 (0). The values for £(0) was obtained as 
1.68 nm, 3.87 nm and 4.33 nm from the plots obtained for T c on , T™ %d and , respectively. 
In order to make an estimate of the lower critical field, H c \, the virgin curves (corresponding 
to the isothermal M-H measurements) were obtained after zero field cooling. The criteria 
of deviation from linearity within virgin M-H plots was used to determine H cl . In Fig. 3, 
we plot the H cl data for both H || (101) and H 1 (101)). A BCS fit (28] yields a H cl (0) 
of 67 Oe and 33 Oe for H || (101) and H _L (101)), respectively. The values of penetration 
depth A(0), Ginzburg-Landau parameter k(0) along with the critical fields are tabulated 
in Table I for FeCro.ccSe. A comparison is made with the superconducting parameters of 
related Fe-based superconductors. Is should be noted here that the BCS fit is not very good 
at temperatures close to T c (0) where the experimental data shows a curvature reminiscent 
of multi-band superconductors. 

In case of type-II superconductors, it is feasible to obtain a practical estimate of the 



critical current density (J c ) via suitable magnetization measurements 29l l30j. This is useful 
when the contact resistance in a typical magneto-transport measurement is such that passing 
large currents through the sample for estimation of J c is impractical without resorting to 



special protocols 



21] for making contacts on the sample. In our case we resort to the 



contactless technique via magnetization measurements for determination of J c (H, T) within 
Bean's model 0, [k]]. Typical five quadrant isothermal magnetization hysteresis (M-H) 
loops recorded at several temperatures between 2 K and 8 K (data at other temperatures not 
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shown here for clarity) are shown in Fig. 4(a) (for H || (101)) and Fig. 4(b) (for H _L (101)). 
The peak in magnetization located at a field value little above nominal zero field in a given 
M-H loop represents the first magnetization peak characteristic, which amounts to (near) 
full penetration of the applied field in the bulk of the sample after zero field cooling. We 
do not observe any prominent feature corresponding to fishtail effect /second magnetization 
peak in our crystals. The absence of a magnetocrystalline anisotropy is evident here too 
(see Fig. 2, where, the values of H C 2for both H || (101) and H _L (101) are presented). 



Making use of the Bean's critical state model formalism 29|, l30j . we extracted critical 

current density values from the isothermal M-H data. Within Bean's model, J c = 20 ,^ M a \ 

00(1-35) 
, with a (= xx mm) and b (= xxmm) being the sample dimensions perpendicular to the 

field direction, and AM is the difference between the magnetization measured during the 

return and the forward legs of the M-H loop. Figure 5 shows the plot of the critical current 

density J C (H) at temperatures between 2K and 10 K for both H || (101) (panel (a)) and 

H 1 (101) (panel (b)). 

Vortex pinning in type-II superconductors can be broadly classified into two categories, 

viz., 5T C pinning (arising because of the spatial fluctuations in the transition temperature, 

T c across the sample) and 51 pinning (caused by the spatial variations in the charge carrier 

mean free path, /). In the case of 5T C pinning, the normalized critical current density, 

-1/2 



J c (t)/J c (o) = (i-t 2 )!/ 6 

where t = T/T c (0) [31 



(1 + 1 2 ) 5 / 6 , while for 51 pinning, J c (t)/J c (0) = (1- t 2 ) 5/2 (l + 1 2 
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321 ]. Following a recipe similar to Das et al 19[ and Thakur et al 



34j, in Fig. 6, we make a comparison of the experimentally obtained J c values (for H 



(101)) with the theoretically expected variations within the 5T C and 51 pinning scenarios. The 
observations point to the dominance of the 51 pinning mechanism in FeCr .o2Se, suggesting 
the occurrence of single vortex pinning by randomly distributed weak pinning centers. 

In the presence of an external magnetic field, the dissipation behavior of a typical type-II 
superconductor is guided by the competition between the pinning force (due to 5\ or, 5T C 
pinning) and the Lorentz force (acting on the flux line because of the transport current). 
The dissipation behavior comprises of two distinct regimes: (a) the pinning force dominated 
flux-creep regime, and (b) the Lorentz force dominated flux-flow regime. In the case of 



high T c superconductor Bi2.2Sr 2 Cao.8Cu20 8 , Palstra et al |35|, |36[ conceived of a dissipation 
behavior guided uniquely by an orientation and field dependent activation energy Uq{H, <p). 
It is of interest to evaluate the extent to which such a formalism is applicable for the Fe- 
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based superconductors. The main panel in Fig. 7 shows the variation of resistivity (p) as 
a function of temperature for a suitably cleaved single crystal of FeCr x Se at T and 9 T 
(|| (101)). The inset panel (a) shows the variation of 4^ with temperature whereas in inset 
panel (b) is plotted the variation of p (close to the superconducting transition) as a function 
of temperature for H || (101) and H _L (101), respectively. In order to obtain the activation 
energy Uq, in Fig. 8, we plot In p versus 1/T at various fields H for both H || (101) (panel 
(a)) and H _L (101) (panel (b)). Linear fits are performed on the low temperature part (large 
1/T behavior) and these are shown by solid lines in both the panels of Fig. 8. The slopes 
obtained from these fits yield the values of the activation energies Uq at various field values 
for both H || (101) and H _L (101). As we know, for thermally activated flux-flow behavior, 
In p(T,H) = lnpo(H) — Uq(H)/T [35]. The plots of Inp versus 1/T at various fields have a 
common intersection point (see panels (a) and (b) of Fig. 8)which should ideally correspond 
to the T c of the superconducting single crystal. In our case this point corresponds to a 
1/T value of ~ 0.078. This corresponds to a temperature of 12.8 K which is indeed very 
close to the T™ ld for our single crystals of FeCro.o2Se. The variation of Uq with H (for both 
H || (101) and H _L (101)) is shown in Fig. 9. We find very similar behavior for both the field 
orientations. Power law fits (~ H~ a ) are performed on the Uq (H) data and two regimes are 
clearly evident with a crossover occurring at a characteristic field H cross (« 38 kOe) for both 
H || (101) and H _L (101). The values for the exponent a in the two field regions for different 



field orientations are provided in Tab 



e 2. It should be noted here that a similar analysis for 



Bi 2 .2Sr 2 Cao.8Cu208 by Palstra et al [35|, [36] also yielded a similar crossover field behavior 



with different values of a on the two sides. However, in their case, very different values for a 
were obtained for the two field orientations (possibly related to the large anisotropy in single 
crystals of B^^S^Cao.gC^Og). Recently Lei et al [373 have found a similar behavior for 
single crystals of /3-FeSe with H cross « 30 kOe (smaller than the crossover field of ~ 38 kOe 
observed in our single crystals). In the case of /3-FeSe, a ~ 0.25 for H < H cross while for 
H > H cross , a « 0.68 (H || (101)). The values of a for the other field orientation (H _L (101)) 
were seen to be within few percent of the values for H || (101). In our single crystals of 
Cr substituted /3-FeSe, we observe a very similar behavior except that the values of a are 
marginally larger. This is understandable considering m 30% enhancement in T c upon 2% 
Cr substitution at the excess Fe site in /3-FeSe. A change in the value of the exponent a 
by about a factor of 3 across H cross with a relatively small value at low fields suggests a 
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crossover from a single-vortex pinning dominated regime to a collective flux creep regime 



The nature of charge carriers in the normal state just prior to the superconducting transi- 
tion is one of the crucial ingredients that is required to physically understand the mechanism 
of superconductivity in any material system. In order to develop an understanding in this 
direction, we perform both Hall effect studies and the thermal transport measurements (in- 
cluding the determination of Seebeck coefficient (S)) as a function of temperature in a single 
crystal piece of FeCr .o2Se. For the Hall measurements, a single crystal piece of FeCr .o2Se 
was mounted on the measurement puck and then repeatedly cleaved using a scotch tape to 
thin down the specimen. Contacts for the Hall measurements were made in the five probe 
geometry using silver epoxy and 30 micron gold wires. The observed Hall coefficient as a 
function of temperature is shown in panel (a) of Fig. 10. From the data it is evident that 
the Hall coefficient is negative below 220 K pointing to a prominence of electron like charge 
carriers in this region. Fig. 10 (c) shows a plot of S(T) at external magnetic field values of 
0T and 5T, respectively. Below the superconducting transition temperature S(T) is zero as 
the charge carriers are involved in the formation of Cooper pairs. The inset panel in Fig. 
10 (c) shows the enlarged view of the S(T) data between 2K to 20 K. It should be noted 
that the S(T) measurements confirms the T c of the sample of 10.5 K as obtained from the 
magnetization and electrical transport studies (c.f., Fig. 1). In the presence of an external 
field of 5T, the shifting of T c to the lower temperatures is also evident from the S(T) data. 
Above T c , S(T) has a positive value and it crosses zero at ~ 65 K. This change of sign is 
not seen in Hall measurements. Above 65 K the S(T) changes sign and become negative and 
attains maximum magnitude of ~ 26 /iV/K at 109 K. The S(T) data above 109 K suggests 
a change in the type of majority charge carriers. The change in the type of charge carrier 
at this temperature is also observed in the Hall measurement. Thus, from the S(T) data, 
it is clearly evident that this material has both types of charge carriers. A modulation in 
the fraction of the two types of charge carriers as a function of temperature is responsible 
for the non-monotonic variations in S(T). This could be crucial information in theoretically 
understanding the nature of superconducting state in this class of materials. 

To summarize, we have grown single crystals of tetragonal phase of FeCr .o2Se (anti-PbO 
type; space group Pi/nmm) where, 2% Cr is substituted in excess at the Fe site. These single 
crystals have a T c which is 30 % larger than the parent /3-FeSe. The basic superconducting 
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properties are corroborated using magnetic susceptibility (both AC and DC) as well as 
electrical transport measurements. An estimate of J c is obtained using isothermal M-H data 
which is also utilized to establish the nature of pinning mechanism in these single crystals. 
A one-to-one comparison between the experimentally obtained temperature variation of J c 
with the corresponding theoretical estimates for 81 pinning and ST C pinning is performed. 
Based on this, the pinning properties within these single crystals are attributed to the 51 
pinning mechanism which corresponds to the occurrence of single vortex pinning by randomly 
distributed weak pinning centers. From low field magnetization measurements and magneto- 
transport studies, fundamental parameters like H c i, H C 2, k, £(0), A(0) are obtained which 
are fundamental to the understanding of magnetic properties of these single crystals. An 
attempt is made to understand the flux-flow behavior in these single crystals within the 
thermally activated flux-flow model. The results points to a crossover from single-vortex 
winning; dominated regime to a collective flux creep regime across a characteristic field H cross 
which is practically independent of the magnetic field orientation direction. Finally, the 



issue of nature of charge carriers in the normal state of these single crystals is addressed. 
For this, we resort to Hall effect, thermopower and thermal conductivity measurements 
which together suggests a modulation in the fraction of the two types of charge carriers with 
temperature within the normal state. This could prove to be a vital input for theoretically 
understanding the nature of the normal state within these class of superconductors and 
possibly the physics of the mechanism of superconductivity in these systems. 
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TABLE I: Comparison of superconducting parameters of the single crystal of FeCro.o2Se with other 
related compounds. 



Compounds 


T c 


H c i (Oe) 


H c2 (kOe) 


£ (nm) 


A (nm) 


K 


anisotropy 


FeCro.o 2 5e (H || (101)) 


10.5 


66 


220 


3.87 


334 


86.3 


Isotropic 


/3-FeSe 


8.5 


75(1) 


180 


4.28 


~ 309 


72.3 


Isotropic 


FeSei-ajTez 


14 


< 100 


650 


2.2 


560 


254 


3.1(2) 


FeTe . 8 S .2 


8.4 




440 


2.7 






1.05(4) 



TABLE II: Power law decay of activation energy (JJq ~ H a ). 



Sample 


Field Orientation 


H C ross (kOe) 


a 


(H < H cross ) 


(H > H cross ) 


FeCr .o2Se 


II (ioi) 


39.3 


0.37±0.03 


0.88±0.04 


1 (101) 


37.3 


0.31±0.02 


0.85±0.03 


Bi2.2Sr 2 Cao.8Cu 2 8+5 [35, 36] 


a, b 


10.0 


0.48±0.04 


0.15±0.02 


_L a, b 


30.0 


0.16±0.02 


0.33±0.05 


/3-FeSe [37] 


II (ioi) 


30.0 


0.25±0.06 


0.68±0.06 


± (101) 


30.0 


0.26±0.02 


0.70±0.09 


SmFeAsOo.85 [39] 


c 


30.0 


w 0.35 


0.88 


± (101) 









11 



5 



10 



15 




FIG. 1: (Color online) The x-ray diffraction pattern with identification of (00/) lines for a cleaved 
piece of single crystal of FeSeo.sTeo.5 (see text for details). 
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FIG. 2: (Color online) The upper critical field, H C 2(T) obtained using different criteria (see text 
for details) for H || (101) (panel (a)) and H _L (101) (panel (b)). 
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T(K) 

FIG. 3: (Color online) The lower critical field, H cl (T) obtained for both H || (101) and H J_ (101). 
The inset panel shows a portion of the virgin isothermal M-H data for both H || (101) and H _L (101) 
obtained at 5K. 
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FIG. 4: (Color online) Isothermal magnetization (M— H) data obtained at various temperatures 
for H || (101) (panel (a)) and H _L (101) (panel (b)) 
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FIG. 5: (Color online) Critical current density (J C (H)) data obtained at various temperatures for 
H || (101) (panel (a)) and H _L (101) (panel (b)). 
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t( = T/T(0)) 

FIG. 6: (Color online) Theoretically expected variation of the normalized critical current density 
(J c ) with reduced temperature within the 51 (solid black line) and <5T C (dotted red line). Also 
plotted are the experimantally obtained data points for (J c ) obtained at typical magnetic field 
strengths of 5 kOe and 10 kOe for H || (101). 
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FIG. 7: (Color online) Variation of resistivity (/?) as a function of temperature for a suitably cleaved 
single crystal of FeCr x Se obtained in an external dc magnetic field of T and 9 T, respectively for 
H || (101). Inset (a): Variation of with temperature. Inset (b): Variation of p as a function of 
temperature for H || (101) and H _L (101) respectively. 
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FIG. 8: (Color online) lnp versus 1/T at various externally applied field values for (a) H || (101) 
and (b) H _L (101). 
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FIG. 9: (Color online) Uq versus H for H | (101) and H _L (101). The power law fits are shown 
by bold lines. 
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FIG. 10: (Color online) (a) Variation of Hall coefficient (Rh) with temperature of a suitably 
cleaved single crystal of FeCr^Se (see text for details) in the presence of an external magnetic field 
of strength 5 T for H || (101). (b) Variation of thermopower (S) with temperature for a single 
crystal of FeCr^Se at T and 5 T, respectively 2dr H || (101). Variation of S with temperature near 
T c is shown as an inset in panel (b). (c) Variation of thermal conductivity (k) with temperature 



